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ABSTRACT 

The  application  is  reported  of  a  United  Sensor  Corpor- 
ation 5-hole  pressure  probe  to  determine  the  flow  field 
upstream  and  downstream  of  the  rotor  of  a  transonic  com- 
pressor operating  to  near  sonic  speeds.  Calibration  of  the 
probe  was  reduced  to  a  set  of  equations  to  allow  off-line 
data  reduction  and  presentation .using  a  Hewlett-Packard 
Model  9830A  Calculator  system.  The  flow  velocities  and 
angles  measured  upstream  of  the  rotor  were  compared  with 
results  of  compressible  flow  calculations.  Blockage  factors 
were  determined  which  related  analytical  approximations 
of  the  flow  incident  at  the  rotor  to  the  measured  flow 
rate.  The  measurements  downstream  of  the  rotor  were  com- 
pared with  those  of  a  second  probe  of  different  design. 
Both  probes  required  corrections  to  their  calibrations 
which  were  established  in  an  axial  annular  flow  of  sim- 
ilar geometric  cross  section.  It  is  shown  that  the  results 
derived  from  the  corrected  calibration  of  the  5-hole  probe 
were  the  more  accurate. 
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I.   INTRODUCTION 

The  study  reported  here  is  an  initial  part  of  a  program 
to  obtain  a  full  understanding  of  the  flow  within  a  tran- 
sonic compressor  stage,  with  the  ultimate  object  of  obtaining 
improved  designs.  Toward  this  goal,  a  transonic  compressor 
of  the  "impulse"  type  was  designed  and  built  by  Dr.  M.H. 
Vavra  in  the  Turbopropulsion  Laboratories  at  the  Naval 
Postgraduate  School.  The  work  was  sponsored  by  Naval  Air 
Systems  Command,  Code  310  through  the  office  of  Mr.  H.J. 
Mueller.  The  so-called  TRANSX  compressor  has  been  operated 
to  6$%   of  design  speed. 

The  inlet  of  the  TRANSX  compressor  is  designed  to  con- 
tract and  to  rapidly  accelerate  the  flow  in  order  to  'gener- 
ate near  uniform  conditions,  with  thin  boundary  layers,  at 
the  rotor  face.  It  is  necessary  to  determine  the  radial  dis- 
tributions of  the  flow  velocity  and  angle  in  order  to  estab- 
lish the  operating  condition  and  to  evaluate  the  performance 
of  the  rotor  blading.  At  or  near  design  speeds,  the  possibil- 
ity that  wakes  in  the  flow  might  induce  vibrations  in  the 
highly  stressed  rotor  blades,  prohibits  the  use  of  probes 
just  upstream  of  the  rotor.  Therefore  it  is  necessary  to  be 
able  to  predict  velocity  and  pitch  angle  distributions  at 
high  speeds  from  data  gathered  at  lower  speeds.  Reference  1 
gives  analytical  approximations  of  the  rotor  face  flow  velo- 
city and  angle  distributions  which  were  obtained  from  com- 
pressible flow  calculations.  The  first  goal  of  the  present 
study  was  to  determine  the  flow  into  the  rotor  by  measurement, 
and  to  compare  the  results  with  the  analytical  predictions. 
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A  United  Sensor  5-hole  cylindrical  pressure  probe  was  selec- 
ted for  its  ability  to  measure  pitch  angle  and  for  ease  of 
mounting.  The  probe  calibration  was  first  established  in  a 
uniform  jet.  Following  a  description  of  the  TRANSX  compress- 
or in  Section  II,  the  investigation  upstream  of  the  rotor  is 
described  in  Section  III  of  this  report. 

Measurement  of  the  flow  downstream  of  the  rotor  is 
necessary  in  order  to  establish  blade  element  performance. 
The  calibrated  United  Sensor  5-hole  probe  and  a  NASA  combin- 
ation probe  [Ref.  2]  were  used  simultaneously  to  conduct 
radial  surveys  of  the  flow  downstream  of  the  rotor.  In  this 
manner  redundant  measurements  of Mach  number  and  flow  angle 
were  obtained  using  two  completely  independent  probes  and 
their  respective  calibrations.  It  was  found  that  the  Mach 
numbers  evaluated  by  the  two  probes  disagreed  by  as  much  as 
10fo   when  the  calibrations,  established  in  the  uniform  jet, 
were  applied  to  the  probe  measurements  made  downstream  of 
the  rotor.  Methods  were  devised  for  correcting  the  calibra- 
tions established  in  the  uniform  jet  for  application  of  the 
probes  in  the  smaller  compressor  annulus.  Following  the 
application  of  corrections,  which  were  greatest  for  the  NASA 
probe,  the  disagreement  in  Mach  number  profiles  was  greatly 
reduced;  in  some  cases  to  less  than  1%,   The  measurements 
carried  out  downstream  of  the  rotor  using  the  5-hole  probe 
are  described  in  Section  IV. 

The  calibration  of  the  5-hole  probe  and  the  method  de- 
vised for  correcting  the  calibration  for  application  in  the 
compressor  annulus  are  described  fully  in  Appendix  A. 
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Inter- related  non-linear  calibration  relations  are  therein 
represented  by  analytic  functions  which  allow  computer  re- 
duction of  probe  survey  measurements.  The  methods  described 
in  Appendix  A  can  be  adapted  to  other  probes  and  to  measure- 
ments with  the  same  probe  in  passages  of  different  geometry. 

II.   THE  TRANSONIC  COMPRESSOR 

A.  DESCRIPTION 

A  schematic  of  the  compressor  test  rig  is  shown  in  Fig. 
1.  The  test  rig  consists  of  an  air-turbine  drive  unit  and 
an  induction  section  containing  a  filter,  a  throttle,  a 
settling  chamber  and  a  flow  measuring  nozzle.  The  air  tur- 
bine is  designed  to  deliver  ^50  H.P.  at  30,000  RPM.  The 
transonic  compressor  is  designed  to  operate  at  30,^60  RPM 
with  a  relative  tipMach  number  of  1»5«  The  compressor  flow 
rate  is  controlled  by  a  rotating  throttle  plate  L^ige  l]« 
The  test  rig,  located  in  an  explosion  proof  cell,  is  con- 
trolled by  an  operator  at  a  remote  console.  A  complete  des- 
cription of  the  laboratory  facilities  ere  given  in  Ref.  3» 
and  the  details- of  the  test  rig  design  is  given  in  Ref.  ^. 
Figure  2  shows  the  compressor  assembly  and  identifies  the 
instrumentation  stations  0  through  3  referred  to  in  the 
following  discussion.  As  indicated  in  view  J  of  Fig.  2, 
peripheral  surveys  were  made  by  rotating  the  outer  casing, 

B.  INSTRUMENTATION 

The  transonic  compressor  is  equipped  with  stationary 
probes  to  provide  routine  performance  data.  The  measurement 
of  torque  was  not  available  at  any  time  during  this 
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investigation,  since  the  honeycomb  flow  straightener  (Fig. 
3)  was  removed  to  allow  higher  flow  rates.  The  fixed  instru- 
mentation is  not  discussed  in  this  report. 

The  results  presented  here  were  obtained  using  a  United 
Sensor  DA-125  three  dimensional  pressure  probe.  The  probe 
calibration  is  reported  in  detail  in  Appendix  A.  The  probe 
mounting  arrangement  is  shown  in  Fig.  3«  Radial  surveys  were 
'taken  at  a  fixed  peripheral  angle  by  manually  positioning 
the  probe  using  a  traverse  unit.  Peripheral  surveys  were 
made  by  manually  rotating  the  outer  casing,  causing  the 
probe  to  move  through  an  arc  at  a  fixed  radius.  All  press- 
ures were  measured  by  a  calibrated  transducer  in  a  scani- 
valve  arrangement.  Predetermined  angular  and  radial  posi- 
tions were  input  into  the  data  logging  system  on  thumb 
wheels. 
C.   DATA  ACQUISITION  AND  REDUCTION 

A  data  logging  system  [Ref .  3]  was  used  to  scan  and  re- 
cord pressure,  temperature  and  position  data.  The  recorded 
output  was  on  punched  paper  tape.  Reduction  of  data  was 
carried  out  using  a  Hewlett-Packard  Model  983OA  Calculator 
system,  shown  schematically  in  Fig.  ty, 

A  program,  designated  TCR  Reduction  (Station  1),  was 
written  to  reduce  data  from  rotor  inlet  flow  measurements. 
Data  were  input  from  the  punched  tapes.  The  program  calcula- 
ted and  printed  the  non-dimensional  radial  probe  position, 
compressor  case  angle,  flow  pitch  and  yaw  angles,  non-dimen- 
sional velocity,  velocity  and  Mach  number.  These  results 
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were  stored  on  digital  magnetic  tape  and  the  radial  distri- 
butions were  graphed  using  the  X-Y  plotter  (Fig.  4). 

A  program  entitled  TCR  Reduction  (Station  2)  was 
written  to  calculate  the  flow  downstream  of  the  rotor  from 
probe  measurements.  The  program  accepted  data  from  magnetic 
tape,  rather  than  punched  paper  tape.  The  paper  tape  data 
were  therefore  transferred  to  magnetic  tape  before  beginning 
data  reduction.  This  procedure  was  selected  because  data 
from  two  different  probes  were  required  to  be  in  the  same 
form  to  allow  distribution  plots  to  be  made  on  a  common 
graph.  When  data  had  been  reduced,  the  calculated  results 
were  stored  on  magnetic  tape  and  the  radial  profiles  were 
produced  on  the  X-Y  plotter. 

III.   FLOW  INTO  THE  ROTOR 
A.   DESCRIPTION  OF  TESTS 

The  purpose  of  these  surveys  was  to  measure  the  flow 
properties  at  the  rotor  inlet  and  to  compare  the  results 
with  analytical  distributions  reported  in  Ref.  1.  Prelimin- 
ary impact  pressure  surveys  were  made  prior  to  this  investi- 
gation, and  as  a  result,  a  bellmouth  inlet  fairing  was  added 
(Fig.  2)  to  improve  inlet  profiles. 

The  5-bole  probe  was  positioned  ahead  of  the  rotor  at 
station  1,  The  flow  profiles  were  measured  in  two  separate 
runs  (numbered  33  and  39).  All  surveys  were  radial,  from 
hub-to-tip.  The  speed  was  set  at  50%   of  design  speed  in  run 
33t  and  varied  from  6$%   to  k0%   of  design  speed  in  run  39» 
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Throttle  settings  corresponded  to  fully  open,  near  surge,  or 
design  point  for  each  speed  of  a  given  survey.  The  case 
angle,  or  peripheral  location  of  the  probe,  was  not  consider- 
ed to  be  significant  for  measurements  at  this  station. 

The  surveys  began  with  the  sensing  holes  0.27  inch 
from  the  rotor  hub.  The  probe  was  then  withdrawn  in  0,1 
inch  intervals  until  the  total  pressure  sensing  port  was 
flush  with  the  inner  surface  of  the  casing.  At  each  point, 
pressure,  temperature  and  position  data  were  recorded  on 
punched  tape.  These  data  were  reduced  to  velocities,  angles 
and  non-dimensional  position. 
B.   RESULTS 

Radial  distributions  of  the  pitch  angle  and  Mach  number 
derived  from  probe  measurements  are  shown  in  Fig.  5  to  Fig. 
11.  The  measured  data  were  reduced  using  the  probe  calibra- 
tion described  in  Appendix  A,  section  A2,  The  radial  distri- 
bution of  pitch  angle,  derived  in  Ref,  1,  is  also  shown  on 
the  plots.  The  data  were  taken  in  two  separate  tests.  Re- 
sults from  run  33  in  Figs.  5-7   show  good  agreement  with  the 
derived  profile  except  at  the  upper  and  lower  20$  of  the  dis- 
tribution. Note  the  wall  effect  on  pitch  angles  for  radial 
positions  greater  than  R=0,85.  Within  the  measurement  accur- 
acy, distribution  of  pitch  angle  is  independent  of  flow  rate 
for  all  surveys  taken  in  a  given  test.  The  results  of  run  39 
(Figs.  9-12)  show  an  overall  reduced  level  of  pitch  angle,  but 
also  shows  no  dependence  upon  throttle  position  or  speed. 
The  reason  for  the  reduced  levels  is  net  known  and  should 
first  be  confirmed  by  further  measurements. 
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In  order  to  plot  the  analytical  Mach  number  distribu- 
tion as  a  comparison  with  the  measurements  plotted  in  Fig,  5 
to  Fig.  11,  the  area  blockage  factor  must  first  be  determin- 
ed. Therefore  the  measured  distribution  of  the  mass  flux  com- 
pared with  the  analytical  distribution  was  chosen  to  repre- 
sent the  same  information,  and  is  presented  in  the  following 
section. 

The  calibration  correction- for  application  in  the  com- 
pressor annulus  does  not  apply  to  measurements  at  the  rotor 
inlet.  However,  this  correction  was  applied  to  the  profile 
in  Fig.  6,  and  the  results  displayed  in  Fig.  8,  The  correc- 
tion applied  near  the  outer  wall  causes  the  profile  to 
better  match  the  analytical  distribution  at  the  wall.  Since 
the  wall  anglp  at  the  inlet  is  the  same  as  that  in  the  cali- 
bration annulus,  the  effect  brought  about  by  applying  the 
correction  supports  its  validity  for  this  region. 
C.   DETERMINATION  OF  BLOCKAGE  FACTOR 

Details  of  blockage  factor  calculations  are  given  in 
Appendix  B.  Reference  1  determined  a  viscous  blockage  factor, 

K„  =0.989*  which  accounted  for  wall  boundary  layer  displace- 

Bl 
ment  thickness.  Figures  13  and  1^  show  distributions  of  the 

non-dimensional  mass  flux,  or  flow  function  obtained  from 
probe  measurements  and  a  comparison  with  the  analytical  dis- 
tribution given  in  Ref.  1  and  Appendix  B.  It  is  shown  in 
Appendix  B  that  by  equating  the  integral  of  the  measured  mass 
flux,  over  the  region  for  which  the  probe  measurements  are 
free  of  wall  effects,  to  the  integral  of  the  analytical  mass 


21 


flux  distribution,  a  blockage  factor  can  be  defined.  The 
blockage  factors  listed  in  Figures  13  and  1*4-  were  obtained 
in  this  manner  and  the  analytical  distributions  are  consis- 
tent with  these  values. 

The  prediction  that  the  mass  flux  increases  with  in- 
creasing radius  is  supported  by  the  measurements,  but  the  mea- 
sured  variation  is  larger  than  anticipated.  This  could  be 
the  result  of  the  inlet  bellmouth  velocity  profil-e  which  has 
yet  to  be  measured. 

The  mass  flux  profiles  are  quite  similar  in  shape,  but 
considerable  variation  occurs  in  the  blockage  factors  de- 
duced from  the  measurements  (.963  to  .98*0.  This  is  possibly 
due  to  the  measurement  of  the  total  mass  flow  rate,  which 
was  observed  to  vary  during  the  surveys.  As  shown  in  Appen- 
dix B,  an  attempt  was  made  to  eliminate  this  effect;  however 
insufficient  data  was  taken  to  make  the  process  statisti- 
cally accurate.  It  is  observed  that  the  smallest  blockage 
factors  were  obtained  at  the  lowest  flow  rates,  for  which 
the  measurements  are  least  accurate.  However,  it  is'  also 
noted  that  at  highly  reduced  throttle  settings,  greater  dis- 
tortion appears  in  the  profiles  near  the  rotor  tip,  suggest- 
ing that  an  increase  in  blockage  does  indeed  occur.  This 
question  may  be  resolved  by  collecting  more  data  with  the 
objective  of  determining  how  blockage  varies  with  RPM  and 
throttling  conditions. 

Figures  15  and  16  show  calculated  flow  distributions 
into  the  rotor  using  the  method  of  Ref.  1  with  a  blockage 
factor  of  0.98.  Note,  in  Fig.  15,  that  choking  of  the  inlet 
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occurs  at  slightly  greater  than  21  lb/sec.  As  blockage  in- 
creases (decreasing  effective  cross  sectional  inlet  area) 
choking  occurs  at  lower  flow  rates.  Figure  16  shows  that  the 
rotor  blade  relative  Mach  number  is  subsonic  everywhere  at 
the  speeds  encountered  thus  far.  At  design  speed  (30,460 
RPM)  S0%   of  the  blading  will  see  supersonic  relative  veloci- 
ties, assuming  a  blockage  factor  of  approximately  O.98. 

IV.   FLOW  DOWNSTREAM  OF  THE  ROTOR 
A.   DESCRIPTION  OF  TESTS 

The  United  Sensor  probe  and  a  NASA  combination  probe 
[Ref.  2]  were  mounted  at  station  2,  separated  peripherally 
by  a  kS     angle.  There  was  a  1.5°  difference  in  the  position 
of  each  probe  with  respect  to  a  stator  blade  leading  edge. 
Peripheral  surveys  were  plotted  such  that  peripheral  angle 
displacement  of  each  probe  was  measured  with  a  stator  blade 
leading  edge  as  reference.  During  radial  surveys  the  two 
probes  were  moved  outwards  and  data  from  both  were  taken 
simultaneously. 

Three  separate  runs,  numbered  43,  44  and  45  are  re- 
ported. The  tests  consisted  of  one  radial  survey  and  one 
peripheral  survey  at  the  passage  centerline  for  each  of  the 
three  runs.  Runs  43  and  44  were  conducted  at  the  same  RPM, 
but  at  different  throttle  settings.  Runs  44  and  45  were  con- 
ducted at  the  same  throttle  setting,  but  at  different 
speeds. 

The  purpose  of  these  surveys  was  to  verify  measurements 
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through  redundancy  and  to  calculate  rotor  performance. 
Agreement  was  obtained  only  after  probe  calibrations  were 
corrected,  as  described  in  Appendix  A.  The  verification  of 
the  measurements  themselves  is  described  in  this  section,  and 
the  procedure  for  calibration,  correction  and  application  to 
the  5-hole  pressure  probe  is  described  in  Appendix  A.  The 
evaluation  of  rotor  blade  element  performance  will  be  re- 
ported separately. 
B,   RESULTS 

Figure  17  shows  a  comparison  of  Mach  numbers  for  run  ^5» 
point  3 i  evaluated  from  5-hole  and  NASA  probe  measurements 
using  calibrations  established  in  a  uniform  jet.  Approxi- 
mately 10%   disagreement  between  the  measured  Mach  numbers 
was  observed,  Since  neither  the  static  pressure  distribu- 
tions calculated  from  the  5-hole  probe  calibration  nor  that 
from  the  NASA  probe  calibration  were  consistent  with  the 
measured  compressor  hub  and  tip  wall  pressures,  the  Mach 
number  distributions  shown  in  Fig.  17  were  not  acceptable. 
It  had  also  been  observed  from  surveys  of  the  uniform  jet, 
described  in  Appendix  A,  that  the  probes  indicated  uniform 
flow  profiles  over  different  regions  of  the  flow.  In  parti- 
cular, the  NASA  probe  failed  to  indicate  uniform  conditions 
until  it  was  traversed  radially  more  than  3  inches  into  the 
flow.  The  effect  of  "immersion"  or  "wall  effect"  on  the 
calibration  of  the  two  probes  was  the  probable  cause  of  the 
observed  disagreement  exhibited  in  Fig.  17.  Since  the  magni- 
tude of  the  immersion  effect,  in  the  uniform  jet,  was  quite 
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different  for  each  probe,  the  approach  taken  was  to  first 
devise  separate  corrections  to  the  calibrations  of  the  two 
probes.  Agreement  between  the  corrected  survey  measurements 
with  both  probes,  and  consistency  of  the  calculated  static 
pressure  distributions  with  the  measurements  at  the  hub  and 
tip  walls  v/ere  sought.  The  corrections  were  based  upon  sur- 
vey measurements  made  in  an  annular  axial  flow  passage  which 
was  similar  in  size  to  the  compressor  annulus  downstream  of 
the  rotor.  The  procedure  for  correcting  the  5-hole  probe 
calibration  is  described  in  Section  A3  of  Appendix  A. 

The  results  shewn  in  Fig.  18  to  Fig.  29  were  obtained 
following  the  application  of  corrections  to  the  measurements 
of  both  probes.  The  upper  sections  of  Figs„  18-23  show  the 
impact  pressure  (to  the  right)  and  indicated  static  press- 
ure (dashed  lines)  measured  by  the  probes.  The  solid  lines 
toward  the  left  are  the  distributions  of  static  pressure 
evaluated  from  the  probe  measurements  after  the  annulus 
corrections  v/ere  applied.  The  lower  sections  of  Fig.  18  to 
Fig.  23  compare  Mach  number  and  flow  angles  derived  from  the 
probe  measurements.  Figures  24  to  29  show  pitch  angle  dis- 
tributions from  the  5-hole  probe  measurements. 

The  results  shown  in  Figs.  18  to  29  were  from  the  three 
tests  mentioned  above.  Figures  18,  19  and  20  show  radial 
surveys  at  $0%   of  design  speed  and  design  incidence  (deter- 
mined from  Fig,  15)  t  50fo   speed  near  surge  and  65%   speed  near 
design  incidence,  respectively.  Figures  21,  22  and  23  show 
peripheral  surveys  of  the  centerline  of  the  channel  at  the 
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same  test  conditions.  Figures  2k   to  29  show  the  resulting 
pitch  angle  distributions  corresponding  to  tests  shown  in 
Figs.  18  to  23. 
C'  DISCUSSION 

1.  Radial  Surveys 

The  radial  surveys  are  shown  in  Figs.  18  to  20.  The 
Mach  number  comparison  between  the  two  probes  is  in  much 
better  agreement  than  the  uncorrected  results  of  Fig.  17. 
The  differences  close  to  the  hub  are  a  result  of  the  static 
pressure,  rather  than  total  pressure  measurements.  It  was 
noted  that  corrected  static  pressure  measurements  for  the 
5-hole  probe  tended  toward  measured  hub  and  tip  wall  press- 
ures, whereas  the  NASA  probe  corrected  static  pressure  tend- 
ed to  values  that  were  too  high.  As  an  example,  for  the  test 
shown  in  Fig.  20,  the  measured  hub  and  tip  wall  pressures 
were  -5.3  and  23.1  inches  of  water  respectively.  Therefore, 
the  Mach  number  calculated  from  the  5-hole  probe  measurements 
is  considered  to  be  more  accurate  than  that  given  by  the  NASA 
probe . 

Yaw  angle  measurements  of  flow  behind  the  rotor  are  in 
agreement  to  better  than  one  degree.  The  reason  for  the  resi- 
dual disagreement  is  not  clear.  It  was  noticed  that  the 
largest  and  most  consistent  disagreement  occurred  at  the 
highest  speed  and  flow  rate  (Fig.  20).  In  Fig.  21  the  flow 
angle  at  the  rotor  tip  is  large;  corresponding  to  a  low  axial 
velocity  component.  This  suggests  a  stalled  condition  at  the 

rotor  tips. 

Total  pressure  measurements  do  not  agree  as  well  as  in 
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the  steady  flow  in  which  the  probes  were  calibrated  (Appendix 
A).  Leaks  were  thought  to  have  been  eliminated  from  the 
scanivalve  and  pressure  transducer  data  system.  There  is 
a  possibility  that  the  total  pressure  differences  were 
caused  by  the  unsteady  nature  of  the  flow.  This  argument  is 
supported  by  Fig.  19.  corresponding  to  a  near  stall  condi- 
tion, for  which  additional  unsteadiness  might  be  expected. 
However,  at  the  higher  frequencies  associated  with  the 
higher  RPM  in  Fig.  20,  the  disagreement  is  greater.  Such 
would  not  be  the  case  if  the  increased  differences  were 
caused  by  the  unsteady  response  of  the  probes. 

2.  Peripheral  Surveys 

Figures  21  to  23  show  the  results  of  peripheral 
surveys  taken  at  the  passage  centerline,  corresponding  to 
test  conditions  shown  in  Figs.  18  to  20.  In  Figs.  21  and  23, 
the  impact  pressure  difference  is  small  and  uniform  with  re- 
spect to  angular  position.  (For  these  tests,  profiles 
hub-to-tip,  were  also  similar), At  the  near  stall  condition 
in  Fig.  22,  a  variation  in  impact  pressure  with  peripheral 
position  is  evident.  The  two  measurements  agree  only  over  a 
small  portion  of  a  blade  space. 

The  peripheral  variation  in  corrected  static  pressure, 
as  indicated  by  both  probes,  was  quite  pronounced.  With  the 
stagnation  pressures  being  relatively  constant,  a  variation 
in  Mach  number  proportional  to  the  probe  static  pressure 
differences  was  measured.  Inspection  of  these  figures  leads 
to  the  conclusion  that  the  stator  blades  ware  affecting  the 
flow  upstream.  The  minimum  Mach  number  and  flow  angle 
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(Fig.  22)  was  measured  in  line  with  the  leading  edge  of  a 
stator  blade,  at  a  case  angle  of  approximately  3630.  That 
angle  is  consistent  with  the  geometry  of  the  flow  into  the 
stator.  It  is  clear  that  peripheral  variations  in  flow  pro- 
perties must  be  taken  into  account  for  specification  of 
rotor  performance  from  probe  measurements  at  station  2. 
This  is  particularly  true  at  reduced  throttle  settings,  as 
in  Fig.  22. 

The  distribution  of  static  pressure  from  the  NASA 
probe  shown  in  Fig.  22  is  similar  in  shape  to  those  from  both 
probes  in  Figs.  21  and  23.  The  more  uniform  distribution  of 
the  5-hole  probe  static  pressure  is  a  marked  departure  from 
the  distributions  of  Figs.  21  and  23.  Since  the  probes  were 

stator  passage  downstream  of  the  5-hole  probe  contained  a 
separation,  or  local  stalled  condition,  whereas  the  passage 
downstream  of  the  NASA  probe  was  unaffected. 
3.  Pitch  Angle  Measurement 

The  pitch  angle  distributions  measured  downstream  of 
the  rotor  are  shown  in  Figs.  24  to  29  and  correspond  to  the 
surveys  of  Figs.  18  to  23.  The  variation  of  pitch  angle  with 
radius  was  nearly  linear  at  design  throttle  settings  (Figs. 
2k   and  26).  The  pitch  angle  is  indicated  to  be  approximately 
-1  close  to  the  outer  wall,  as  might  be  expected  due  to 
boundary  layer  effects.  These  measurements,  which  rely  upon 
the  annulus  corrections  described  in  Appendix  A,  are  remark- 
ably consistent  with  the  expected  distribution.  At  the  point 
where  probe  measurements  were  taken,  the  hub  wall  angle  was 
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+12°.  Extrapolation  of  the  angle  distribution  to  the  wall, 
in  Figs,  2k   and  26, is  consistent  with  this  value.  The  re- 
sults in  Fig.  25,  for  the  near  surge  condition,  shows  a  re- 
gion of  nearly  constant  pitch  angle  near  the  inner  wall  and 
a  tendency  to  increase  rapidly  toward  the  12°  hub  wall  angle. 
This  is  consistent  with  rake  measurements  of  impact  pressure 
taken  at  the  stator  exit,  which  indicated  a  region  of  sep- 
arated flow  at  the  hub  wall.  This  would  have  the  effect  of 
creating  a  large  pitch  angle  gradient  near  the  hub.  The  lar- 
ger region  of  negative  pitch  angle  near  the  tip  would  be  con- 
sistent with  an  increase  in  wall  boundary  layer  thickness. 
This  region  also  exhibits  a  large  pitch  angle  gradient  when 
compared  to  the  center  $0%   of  the  distribution. 

Peripheral  variations  in  pitch  angle,  shown  in  Figs. 
27  to  29  indicate  little  change  at  or  near  design  cond- 
itions. However,  near  surge (Fig.  28)  the  flow  angle  incr- 
eased on  one  side  of  the  stator  blade  and  decreased  on  the 
other  side.  Though  it  appears  to  be  a  significant  effect, 
the  variation  is  no  more  than  one  degree  either  side  of 
the  mean  line  and  would  not  be  significant  in  rotor  per- 
formance calculations. 
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V.   CONCLUSION 

.  In  this  study  a  5-hole  pressure  probe  was  used  success- 
fully to  determine  the  distribution  of  velocity  and  flow 
angle  upstream  of  a  compressor  rotor.  Blockage  factors  were 
determined  from  measurements,  which  allow  the  use  of  analyti- 
cal expressions  to  relate  the  measured  flow  rate  to  the  in- 
cident flow  properties  at  the  rntor. 

The  same  probe  was  applied  downstream  of  the  rotor. 
The  measurements  were  compared  to  those  of  a  NASA  probe, 
which  did  not  provide  a  measurement  of  pitch  angle.  Part  of 
the  differences  in  the  Mach  number  distributions,  seen  near 
the  hub,  was  probably  for  this  reason.  The  5-hole  probe 
gave  more  accurate  Mach  number  distributions  mainly  because 
the  static  pressure  corrections  for  passage  effects  were 
smaller  than  those  for  the  NASA  probe. 

The  differences  in  the  radial  distributions  of  impact 
pressure  were  partially  the  effect  of  the  difference  in  the 
location  of  the  probes  with  respect  to  the  nearest  stator 
blade  leading  edge.  Although  that  difference  was  only  1.5  » 
at  decreased  throttle  settings  the  peripheral  variations  in 
the  indicated  flow  profiles  became  significant.  In  future 
tests  radial  surveys  should  be  conducted  with  both  probes  in 
the  flow  simultaneously,  at  fixed  speed  and  throttle  setting; 
then  the  probes  should  be  interchanged  and  the  surveys  re- 
peated without  changing  conditions.  The  resulting  plots  would 
differ  only  as  they  are  affected  by  the  slight  peripheral 
displacement  of  the  probes.  Peripheral  survey  variations  must 
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be  accounted  for  if  probe  measurements  at  station  2  are  to 
be  used  for  derivation  of  rotor  performance. 

The  effects  of  unsteady  response  on  probe  measurements 
should  also  be  examined.  It  is  reported  that  the  average  in- 
dicated total  pressure  in  unsteady  flow  can  be  greater  than 
the  time  average  of  the  total  pressure  [Ref.  6],  This  implies 
that  two  different  probes  in  the  same  steady  flow  could 
measure  identical  total  pressures,  and  when  placed  in  un- 
steady flow  might  disagree  as  a  function  of  speed  and  throt- 
tle setting.  While  using  Reference  6,  this  effect  was  calcu- 
lated to  be  small  for  the  probes  and  conditions  in  the  present 
study.  The  differences  in  the  measurements  suggest  that  this 
conclusion  should  be  re-examined. 

The  correction  devised  to  compensate  for  wall  effect  on 
the  5~hole  probe  pitch  angle  measurements  v/as  very  effective 
judging  from  the  plotted  pitch  angle  distributions  at  station 
2.  The  overall  effect  of  correcting  the  5-hole  probe  calibra- 
tion for  application  to  the  specific  flow  geometry  of  the 
transonic  compressor,  points  out  the  necessity  of  careful 
calibration  of  pressure  probes  used  for  making  measurements 
of  internal  flows. 
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Figure  3,   End  View  of  Compressor  Assembly 
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Appendix  A:   Determination  of  Velocity  Vector  From 
5 - Hole  Pressure  Probe  Measurements 

A.l   Method  of  Approach 

The  limited  axial  space  available  for  instrumentation  in  most  turbo- 
machinery  flow  passages  severely  limits  the  types  and  sizes  of  pressure 
probes  which  can  be  used.   The  United  Sensor  DA-125,  5-hole  pressure  probe 
is  designed  for  this  application.  The  probe,  when  calibrated,  will 
measure  flow  velocity,  yaw  angle  and  pitch  angle.  A  major  drawback  of 
the  probe  has  been  the  tedious  labor  required  to  extract  velocity  and  an- 
gle data  from  non-linear,  and  inter-related  calibration  curves.  The 
purpose  of  this  appendix  is  to  document  a  method  for  calibrating  the 
United  Sensor  5-hole  probe  such  that,  in  application,  data  reduction  can 
be  performed  by  computer  at  the  time  of  acquisition. 

Results  are  given  here  for  a  particular  probe  which  was  used,  to 
measure  profiles  in  the  TRANSX.  compressor.   It  was  found  to  be  necessary 
to  correct  the  measurements  made  in  the  compressor  annulus  for  effects 
altributable  to  the  passage  size  and  shape.   The  method  is  therefore  in 
two  parts.  First,  the  probe  was  calibrated  with  reference  to  a  Prandtl 
probe  in  a  jet  contained  within  a  larger  pipe.     .   Analytical  approxima- 
tions were  found  which  were  represented  as  a  function  of  the  radial  position 
in  the  annulus.   These  steps  are  described  in  section  A2  and  section  A3 
respectively,  where  methods  of  applying  the  calibration  with  and  without 
the  passage  correction  are  also  given. 
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A. 2   Calibration  and  Application  Free  of  Wall  Effects 
A. 2.1  Probe  Description  and  Analysis 

The  5-hole  probe  is  shown  in  Fig.  Al.   Characteristics  of  this  and 
similar  probes  are  given  in  the  manufacturer's  catalogue  (8  )•   Briefly, 
the  probe  is  rotated  in  the  flow  until  the  pressure  indicated  by  ports  2 
and  3  (p-,  &  p,)  are  equal.   The  yaw  angle  (ry)  is  then  read  on  a  vernier 
scale  attached  to  the  top  of  the  probe.  (Fig.  A2).   The  average  of  pressures 
p  and  p   (poo>  obtained  pneumatically  through  a  mixing  value)  is  referred  to 
as  the  "indicated  static"  pressure  of  the  probe.  Port  1  (p,  )  measures  "in- 
dicated total"  pressure  and  the  difference  between  the  pressures  p.  and  p 
is  a  measure  of  the  pitch  angle.   However,  the  relationship  of  P9^,  and 
(Pi  -  P  )  to  actual  static  pressure,  total  pressure  and  pitch  angle  must 
be  established  by  calibration  in  a  known  flow.  Indicated  static  pressure  is 
known  to  be  more  subject  to  errors  than  total  pressure,  particularly  when 
the  sensing  holes  are  within  about  ten  probe  diameters  of  a  boundary. 

The  calibration  of  the  probe  is  expected  to  depend  on  flow  Mach  number. 
It  is  more  convenient  howei^er  to  non-dimensionalize  the  velocity  using  the 


"limiting"  velocity,  V  =  /2C  T   (for  a  perfect  gas),  when  the  stagnation 
temperature  is  a  constant  in  the  flow.   The  use  of  limiting  velocity  is 
discussed  in  Appendix  A  of  reference  7  .  The  non-dimensional  velocity  is 
given  by 

X  =  V/Vt  A(l) 

where  V  =  109.62  /T  ft/sec.  for  air  when  the  stagnation  temperature,  T  , 
is  in°R.  As  shown  in  reference  7  >  the  relationships  between  the  static 
and  stagnation  gas  properties  become 
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T  2 

tVi-x 


■2-  =    (1   -  X2)  Y_1  A(2) 

Pt 


Y 


t 

The  relationship  between  Mach  number  and  X   is   given  "by 

?        ? 

x  ■  rr^  A(3) 

The  calibration  of  the  5-hole  probe  can  be  represented  in  the  form  of 
two  functions  of  Mach  number  (or  X)  and  pitch  angle,  0. 

If  p   and  p  are  the  true  stagnation  and  static  pressures,  then  define 

pl  "  p?^ 
F  (X,  0)   =  -± £i  A(k) 

a   '  ■   Pt  -  p 

and 

F0  (x>  «  "  57^  A(5) 

to  be  the  required  calibration  functions,  where  it  is  understood  that  the 

probe  pressures  are  taken  when  p  =  p 

<-  .5  • 

The  value  of  X  is  obtained  from  true  static  and  stagnation  pressures 

using  Eq.  A(2) : 

where  p  =  atmospheric  pressure,  and  the  pressure  differences  can  be 
measured  by  manometers. 
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To  calibrate,  the  probe  must  be  placed  in  a  flow  at  a  known  pitch  angle 

and  the  true  static  and  impact  pressures  measured  by  a  reference  probe. 

From  the  measured  pressures,  the  value  of  X  is  calculated  from  Eq.  A (6)  and 

the  values  of  F  and  F^  from  Eq.  A(^)  and  Eq.  A(5).   However,  in  applying 
9       y 

the  probe  in  an  unknown  flow,  expressions  relating  only  5-hole  probe 
measurements  to  the  unknowns  X  and  0  are  required.  Rearranging  Eq.  A(U) 
and  using  Eq.  A(2), 


^— ^  =  (%   F   (X,  0)  [1  -  J2-]  =  (%   F  (X,  0)  [1  -  (1-X2)  Y"1]    A(7) 
*1       yl       e  yt  pl   R 


Define 


■0  W   -'e 

and,  similarly 


Fi  (X,0)  =  Fa  (X,0)  [1  -  (1-X2)  Y_1]  A(8) 


F^  (X,0)  =  F^  (X,  0)  [1  -  (1-X2)  Y-1]  A(9) 

The  two  functions  F'  and  Fi  can  be  tabulated  as  functions  of  X  and  0  using 

9      0  V 

the  values  obtained  in  the  calibration  for  F  and  F  /  at  each  value  of  X. 

9      0 

If,  in  addition,  it  can  be  shown  that  p  /   =  1.0  to  within  the  desired 
accuracy  over  the  range  of  X  and  0  of  interest,  then,  from  Eq.  A (7) 

F'  (X,  0)  =   £  1  '     23J  A(10) 

9  (Pl  *  PA)  +  PA 


and  similarly 


Fi    (X,   0)   =  t VvT A(ll) 

0  r>         (Pj,    -  PA)    +  PA 


6^ 


In  using  the  calibration  expressed  by  Eq.  A(lO)  and  Eq.  A(ll),  it  is 

necessary  to  interpolate  in  the  tables  established  for  F'  and  Fi  to  find 

values  of  X  and  0  for  which  Eq.  A(10)  and  Eq.  A(ll)  are  simultaneously 

satisfied  by  the  measurements.   Alternatively,  analytical  approximations 

can  be  derived  for  F'  and  Fl   which  can  be  solved  by  iteration.   The  analy- 

9      0 

tical  approximation  technique  was  chosen  here. 

A. 2. 2   Calibration  Tests  in  a  Jet 

The  5-hole  probe  was  calibrated  as  shown  in  Fig.  A3-  An  A.S.M.E.  flow 
nozzle,  7  inches  in  diameter,  was  flanged  between  two  36  inch  lengths  of 
'10  inch  diameter  steel  pipe.   (The  internal  contour  of  the  nozzle  was  chrome 
plated  and  polished  to  a  maximum  roughness  of  32  micro  inches).  Air  was 
supplied  to  the  10"  pipe  at  controlled  flow  rates  from  an  industrial  com- 
pressor (  3)  "via  a  h   ft.  diameter  settling  chamber  and  entrance  bellmouth. 
A  12"  long  bundle  of  §•"  diameter  tubing  removed  any  swirl  remaining  in  the 
flow  from  the  plenum.   The  flow  from  the  pipe  downstream  of  the  flow  nozzle 
exhausted  into  the  laboratory. 

A  Prandtl  probe  was  used  to  obtain  the  reference  measurements  of 
static  and  total  pressures.   The  5-hole  and  Prandtl  probe  tips  were  posi- 
tioned on  opposite  sides  of  the  flow  centerline  at  a  radius  of  one  inch. 
The  measuring  station  was  2  inches  downstream  of  the  nozzle  exit  face.   The 
5-hole  probe  was  mounted  through  an  aluminum  block  as  shown  in  Fig.  A*+. 
The  block  allowed  the  pitch  angle  to  be  set  in  5°  increments  from  -  10°  to 
+  10°  and  in  10°  increments  from  -  ^0°  to  +  U0°  with  no  change  in  the  tip 
location.   All  pressure  differences  were  measured  using  U-tube  water  mano- 
meters. The  stagnation  temperature  was  measured  using  an  iron-constantarn 
thermocouple  probe  in  the  flow  upstream  of  the  flow  nozzle. 
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The  flow  at  the  measuring  plane  w*as  first  surveyed  using  fche  Prandtl 
probe  with  the  5-hole  probe  removed.   The  result  in  Fig.  A5  showed  that  the 
velocity  (X)  was  uniform  over  a  V  diam  core,  and  therefore  that  the  flow 
at  the  reference  probe  tip  and  the  5-hole  probe  tip  were  the  same  to  within 
measurement  accuracy. 

The  procedure  for  the  calibration  tests  was  the  following:  The  flow 
rate  from  the  compressor  was  set  and  allowed  to  stabilize.  The  5-hole  probe 
was  then  inserted  in  turn  at  each  desired  value  of  the  pitch  angle.   Pressure 
differences  were  recorded  by  hand  from  the  manometer  for  the  Prandtl  and 
■5-hole  probes.   The  thermocouple  probe  reading  and  barometric  pressure  were 
also  recorded. 

To  cover  the  range  required  in  the  probe  application,  0  was  varied  from 
-5°  to  +  k0°   in  7  steps.   Nine  values  of  X  were  set  giving  dynamic  pressures 

of  reading  each  manometer  column,  the  pressure  differences  are  measured  to 

±  0.05  inches  of  water.   However,  unsteadiness  in  the  air  supply  over  the 

recording  time  period  can  be  expected  to  reduce  to  overall  accuracy  to  +0.1 

inches  water. 

A.2.3  Analysis  of  the  Results 

The  recorded  data  is  given  in  Table  Al.  The  ratio  (p,/p, )  was  calculated 

to  be  within  0.2%  of  unity  over  the  complete  range  of  X  and  0  measured  and 

hence  was  set  equal  to  unity  in  representing  and  applying  the  calibration 

results.   Using  the  data  in  Table  Al,  values  of  F'  and  FL   were  calculated 

e     0 

using  Eq.  A (8)  and  Eq.  A(9),  where  F  and  F-  were  obtained  from  Eq.  A(k)   and 
Eq.  A(5)  ani  X  was  calculated  using  Eq.  A(6).   The  remaining  task  was  to 
find  analytical  approximations  for  the  values  of  F'  and  F\   as  functions  of 

e     0 

X  and  0,  and  to  devise  an  iterative  computer  solution. 
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F'  and  F*,  were  plotted  as  functions  of  0  for  each  value  of  X.   The 
results  are  shown  in  Fig.  A6  and  Fig.  A7.   Using  polynomial  regression  on 
the  Hewlett-Packard  9330A  calculator,  F'  at  constant  X  was  satisfactorily 

e 

fit  by  a  Uth  order  polynomial  in  0.  F\   at  constant  X  was  satisfactorily 
fit  by  a  2nd  order  polynomial  in  0.   This  representation  however  was  found 
to  be  unsuitable  in  solving  iteratively  for  X  and  0  with  F'  and  FL     given. 
Instead,  expressions  were  found  for 

F0  ■  F0  W  'x  =  const 


as  in  Fig.  A7,  and 

X  ■  X  (V  10  -  const 

Fig.  A8  shows  the  results  plotted  in  the  latter  form.   Quadratic  func- 
tions were  found  to  represent  the  data  satisfactorily. 

Since  data  was  taken  for  9  values  of  X  and  7  values  of  0,  the  two  sets 
of  functions  were  expressed  as 

F^   (0)  =  a.(X)  +  b.(X)  0  +  c.(X).  02  A(l2) 

where  i  =  1,  2  9  and  a.,  b.  and  c.  are  coefficients  which  depend  on  the 

value  of  X ,  and 

X.(F')  =  f  (0)  +  g  (0).  [?']  +  h ,(0).[F']2  A(13) 

j  y    j      j     9    j     y 

where  j  =  1,  2 7  and  f-;,  g.  and  h.  are  coefficients  which  depend  on 

J    J      J 

the  value  of  0. 

Next,  the  coefficients  a.,  b.  and  c.  from  the  9  equations  represented 
by  Eq.  A (12)  were  approximated  by  a  polynomial  expansion  in  X,  and  the 
coefficients  f .,  g.  and  h.  from  the  7  equations  represented  by  Eq.  A(13)  by 


6? 


a  polynomial  expression  in  0.   Thus  coefficients  were  found  for  the  expressions 

a.  (X)  -  A,  +  A0X.  +  A.  X2  +  

l  v  '  1    2  1    3   l 

b.  (X)  =  B1  +  B2X.  +  B3  X2±   + A(lUa) 

c.  (x)  =  c±  +  eye.  +  c3x2  +  

f .  (0)  =  Fx  +  F20 .  +  F302  + 

g.  (0)  =  G1  +  G20^  +  G302  + A(lUb) 

h.  (0)  =  Hx  +  H20  .  +  H302  + ■- 

which  gave  the  best  fit  to  the  9  values  of  a.,  b.,  c.  and  to  the  7  values  of 

ill 

f .,  g.  and  h..   This  results  in  a  single  equation  (Eq.  A (12))  for  Fl    (0)  in 

J     J         J  r 

which  the  coefficients  are  analytic  functions  of  X,  and  a  single  equation 

(Eq.  A(13))  for  X(F')  in  which  the  coefficients  are  analytic  functions  of  0. 

9 

The  coefficients  in  Eq.  A(lU)  are  given  in  Table  A3. 

In  application,  the  probe  measurements  give  values  to  F'  and  FL   through 

Eq.  A(10)  and  Eq.  A(ll).  Equation  A(l2)  and  Eq.  A(l3)  must  be  solved  for  the 

unknown  0  and  unknown  X.   An  initial  estimate  of  the  pitch  angle  was  devised 
in  the  following  way.  A  new  function  of  the  probe  pressure  differences  was 
defined  such  that 

*       e      yi       23 

When  0  was  plotted  as  a  function  of  F1  for  each  value  of  X.  it  was 

X 

observed  that  the  curves  were  nearly  linear  with  similar  slopes.   Therefore 

the  average  of  the  9  values  of  F'  at  each  value  of  0  was  calculated  and  a 

A. 

third  order  polynomial  was  found  to  fit  the  7  average  values  when  0  was  plotted 
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as  a  function  the  average  F'  (Fig.  A9)  •   This  gave 

A. 

0   (F1)  =  J  +  K  (F')  +  L  (F')2  +  M  (F*)3  A(l6) 

-*-    A.  A.  A.  A. 

where  0  is  the  desired  initial  estimate  of  the  pitch  angle  and  F'  is 

-L  A. 

obtained  from  the  probe  measurements  using  Eq.  A(15).   The  coefficients  J, 

K,  L  &  M  are  given  in  Table  Ak. 

A2.U  Application 

When  applying  the  5 -hole  probe  to  measure  an  unknown  flow,  the  pressure 

differences  read  when  p  and  p  are  balanced  are  reduced  to  the  three  factors 

F',  Fi  and  F'  using  Eq.  A(10),  Eq.  A(ll)  and  Eq.  A(15).  X  and  0  are  evaluated 
■  0  tP  \ 

using  the  calibration  relations  derived  in  the  preceeding  section.  The  yaw 
angle  is  read  directly  from  the  probe  scale.  The  velocity  is  obtained  from 
the  product  of  X  and  V  which  requires  the  input  of  the  stagnation  temperature. 

The  reduction  to  velocity  and  pitch  angle  was  programmed  for  the  Hewlett- 
Packard  Model  9830A  calculator.  A  flow  chart  of  the  program  is  shown  in  Fig. 
A10.   The  program  first  calculates  the  initial  estimate  of  0  (=0.. )  using 
Eq.  A(l6).   The  initial  estimate  of  X  (=X^ )  is  calculated  from  Eq.  A(l3)  using 
0n  in  Eq.  A(l^b)  to  evaluate  the  coefficients  f ,  g  and  h.  The  coefficients 
a,  b  and  c  in  Eq.  A(12)  are  then  calculated  using  X,  in  Eq.  A(lUa).   In 
general  the  value  of  Fl   which  could  be  calculated  using  these  coefficients  in 
Eq.  A (12)  would  not  be  equal  to  the  value  calculated  from  the  input  data.  An 
iterative  procedure  is  begun  using  0  and  X  as  initial  values. 

In  the  adopted  procedure,  Eq.  A (12)  is  viewed  as  an  equation  to  be  solved 
for  an  unknown  0  in  terms  of  a  known  FL  and  known  X.  Equation  A (13)  is  treated 
as  an  equation  for  X  in  terms  of  a  known  0.  Equation  A (12)  is  solved  by  fixed 
point  iteration  (9)  beginning  with  the  initial  values  0  and  X  .   In  this 
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method,  the  (i+l)   value  of  0  in  the  iteration  is  given  in  terms  of  the  i 
value  by  rewriting  Eq,.  A (12)  as 

0i+1  =  [F^  -  a  -  c  0^1/b  A(17) 

For  a  fixed  value  of  X(=X.),  which  fixes  the  coefficients  a,  b,  and  c, 
successive  values  of  0  are  calculated  using  Eq.  A(17)  until  |0.,,  -0-  ]<£/ 

is  satisfied.   g,  is  a  small  quantity.   A  new  value  of  X(=X.  ,  )  is  then 
p  i+l 

calculated  using  0=0    =0.  in  Eq.  A(13)  and  the  process  repeated  until 
|X.  n  -  X. |  <  g  ,  where  g  is  a  small  quantity. 

The  steps  in  the  solution  may  be  summarized  as  follows: 


Calculate  0  from  Eq.  A(l6) 

Calculate  X  from  Eg.  A(13) 

Calculate  0    from  Eq.  A (17) 

If  |0.  n  -  0.1  <e  /  do  step  5),  otherwise  so  to  step  3) 
i  •  j-    i    p 

Calculate  X.  ,  ffom  Eq.  A (13)  using  0=0  ,  =  0 
i+l  i+l    i 

If  IX    -  X. |  <  e,  stop,  otherwise  go  to  step  3)« 


70 


A. 3  Corrections  to  the  Calibration  for  Application  in  an  Annulus 
A3-1  Introduction 

The  calibration  described  above  was  carried  out  where  the  presence  of 
flow  boundaries  would  not  be  expected  to  influence  I; he  probe  readings.   A 
change-:  in  the  calibration  of  the  probe  is  expected  when  measurements  are 
made  within  about  ten  probe  diameters  of  a  boundary,  and  when  the  probe  size 
is  significant  conipared  to  the  flow  passage  area.  A  precise  correction  which 
would  allow  the  use  of  the  probe  in  any  flow  field  is  a  difficult  task  and 
not  attempted  here.   Instead,  a  calibration  of  the  probe  in  an  axial  annular 
flow  of  similar  geometry  to  the  compressor  annulus  was  carried  out,  and 
corrections  to  the  "wall-free"  calibration  of  section  A2  were  derived  which 
were  represented  as  a  function  of  radial  position.  The  test  of  the  success 
of  this  proceedure  was  to  compare  measurements  made  in  the  compressor  using 
two  different  types  of  probes.  Agreement  in  the  velocity  profiles  derived 
from  the  two  probes  would  be  expected  if  the  corrections  applied  to  each 
calibration  were  valid.  Also,  the  distribution  of  true  static  pressure  de- 
rived from  the  probe  calibrations  should  be  consistent  with  wall  tap  measure- 
ments on  the  tip  and  hub  surfaces.   In  the  following  sections,  the  calibration 
tests  are  described  and  the  derivation  of  corrections  to  the  "wall-free" 
calibration  is  presented.  The  application  of  the  revised  calibration  to 
compressor  measurements  is  then  given. 

A 3« 2  Calibration  Tests  in  an  Annulus 

An  annular  axial  flow  was  generated  as  shown  in  Fig.  All.  An  axisym 
metric  body  with  a  maximum  diameter  of  6"over  28" of  its  length  was  supported 
in  the  center  of  a  10"  pipe.   At  the  upstream  end,  four  thin  airfoil  struts 
were  positioned  ahead  of  the  contoured  nose  section  which  generated  a  smooth 
contraction  of  tne  Plow  passage,  to  a  2"  annulus.   Heavy  struts  supported  the 
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centerbody  at  the  outlet  f  Lange  where  the  flow  exhausted  into  the  laboratory. 
A  view  of  the  apparatus  is  shown  in  Fig.  A.12.   Six  equally  spaced  wall 
pressure  taps  were  provided  in  the  10"  pipe,  ab  each  of  3  axial  stations. 
Corresponding  pressure  taps  were  provided  on  the  centerbody  surface.   Pro- 
vision for  radial  probe  surveys  was  made  at  the  same  three  axial  stations. 

It  was  established  in  preliminary  surveys  that  the  wall  pressures  at 
station  1  in  Fig.  All  were  uniform  around  the  periphery  to  +  0.3Mo  °?  ^ae 
average  dynamic  head.   The  outer  (tip)  wall  pressures  were  consistently 
lower  than  inner  (hub)  wall  pressures  by  0.^+5%  of  the  dynamic  head.   The 
flow  was  measured  to  be  axial  to  within  the  sensitivity  of  the  5-hole 
probe  yaw  angle  measurement  (+  .1°).   Probe  calibration  was  therefore  carried 
out  at  station  1.   Static  pressure  was  assumed  to  vary  linearly  between  hub 
and  tip  values,  and  the  pitch  angle  was  taken  to  be  zero. 

Radial  surveys  were  made  from  hub  to  tip  at  fixed  flow  rates  using 
the  5-hole  probe,  a  Kiel  probe  and  a  pitot-static  probe  in  succession. 
Readings  were  taken  at  the  same  radial  locations  with  each  probe.   Intervals 
of  0.2",  decreasing  to  0.05"  near  the  walls,  were  used.   Dynamic  pressure  was 
varied  from  27  to  73  inches  of  water  in  7-8  inch  increments.  The  data  for 
seven  surveys  is  given  in  Table  A2 . 

A3. 3  Analysis  of  the  Results 

An  example  of  the  distribution  of  pressures  measured  by  the  three  probes 
is  given  in  Fig.  A13.   It  was  observed  that  the  probes  gave  the  same  values 
of  impact  pressure  at  corresponding  points  to  within  the  accuracy  of  the 
measurements  except  close  to  the  outer  wall  where  the  velocity  gradients  were 
large.  The  static  pressure  indicated  by  the  pitot  static  probe  was  observed 
to  vary  very  little  across  the  annulus,  and  was  in  agreement  with  a  linear 
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variation  of  the  true  static  pressure  between  hub  and  tip  measurements 

at  the  wall.  Also,  the  wall  pressures  at  hub  and  tip  away  from  the  traverse 

station  were  not  measurably  affected  by  the  prove  movement.   Consequently,  the 
5-hole  probe  indicated  total  presure  (Pj)  and  a  static  pressure  obtained  by  a 

linear  interpolation  between  measurements  at  the  walls  were  used  to  calculate 
the  actual  velocity  (X)  at  the  probe  location.   By  applying  the  methods  of 
section  A2,  the  5-hole  probe  pressure  measurements  were  reduced  to  "wall-free" 
values  of  velocity  and  pitch  angle,  (X  and  0),  which  in  general  were  different 
from  the  actual  values,  (X  and  0)  at  the  probe  location. 

A  velocity  correction  factor,  K  ,  was  defined  and  examined  as  a  function 
of  position  in  the  annulus  and  the  indicated  velocity,  thus: 

t      (H,X)  =  ^-X-  A(18) 

X  X 

where  H  =  y/h,   y  is  the   distance  outwards   from  the  hub  wall  and  h    is  the 

annulus  depth.     At  each  value,   H. ,   of  the  probe  position  the  dependence   of 

K     on  X  was   approximated  by  a  polynomial 

A 

K        =  q.    +r.X+s.X2+  A  (19) 

x.         ^1  i  i 

l 

where  i  =  1,  2,  —  12,  and  H.  =  0.125,  0.15,  —  0.975-   The  coefficients 
in  Eq.  A(19)  are  listed  in  Table  A5. 

The  pitch  angle  at  each  probe  position,  H. ,  was  found  to  be  independent 
of  X.  A  pitch  correction  factor,  K,  was  defined  which  was  then  a  function 
of  position  only:  thus 

h  (>■)  -  ('*3:g*p)-  *<«,) 

where  p  is  a  constant  angle  introduced  to  avoid  division  by  zero.   Since 
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0  =  0  in  the  calibration  test,  Eq.  A(20)  with  p  =  l80°  becomes 

K^(H)  =  -  |o  A(21) 

The  values  of  K^  obtained  for  the  tests  in  the  annulus  using  Eq.  A (21)  were 
approximated  by  a  6   order  polynomial, 

K,  (H)  =  z±   +  z2  H  +  z3  H2+  zk   H3+z5H1++  z6H5+  z  H6         A(22) 

to  .995  root  mean  square.   The  coefficients  in  Eq.  A(22)  are  listed  in 
Table  A6. 

A3.*-<-  Application 
It  is  assutied  that  the  corrections  required  to  the  calibration  free  of 
wall  effects  when  the  5-hole  probe  is  used  to  determine  the  flow  field  in  a 
compressor  annulus,  are  those  established  for  an  axial  flow  in  the  annulus 
of  similar  geometry.   The  probe  readings  are  first  reduced  to  values  of  X 
and  $  using  the  methods  given  in  Section  A2.   Then,  from  Eq.  A(l8)  and 
Eq.  A (20)  with  p  =  l80°,  the  corrected  velocity  and  pitch  angle  are  given  by 

x  =  (r- V)-  X  A(23) 

1  -  Kx 

and 

0  =  (r-^-r-).  *  +  18°  Cr-V)         A<2U) 
i  -  fy  1  ~  *£ 

At  any  probe  position,  H,  the  value  of  K,   in  Eq.  A(2*+)  is  given  by 
Eq.  A (22)  for  which  the  coefficients  are  given  in  Table  6. 

The  value  of  K  in  Eq.  A (23)  must  be  obtained  by  interpolation.   The 
uncorrected  velocity  X  is  known  at  the  probe  position  H;  however,  the  correc- 
tion factor  K  can  be  evaluated  using  Eq.  A(l9)  only  at  discreet  locations 
H. .   The  method  adopted  is  a  form  of  quadratic  interpolation  for  uneven 
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intervals  (ID).  First  the  value  of  i  is  found  such  that  H  is  located  in  the 

interval  "between  H  and  H    as  shown  in  Fig.  A(lU).   The  values  of  K    , 
•_  i     1+1  xi_1 

K   ,  K     and  K     are  calculated  from  Eq.  A(l9)  using  the  coefficients 

xi   xi+l      xi+2 
given  in  Table  A5.   The  coefficients  of  quadratics  for  K  passing  through 

each  of  the  two  sets  of  three  points  H.  , -  H. ,  H. ,,  and  H. ,  H. ,,,  H.  _ 

l-l   i   i+l      i   i+l'   i+2 

are  calculated  and  each  quadratic  is  used  to  calculate  a  value  for  K  at  H. 
The  arithmetic  mean  of  the  two  values  is  taken  for  the  interpolation.   In 
the  first  or  the  last  interval,  or  to  extrapolate  outside  the  given  range, 
only  one  quadratic  is  used. 

The  data  reduction  with  corrections  was  programmed  for  the  Hewlett- 
Packard  Model  9830A  calculator.   The  flow  chart  is  given  in  Fig.  A(l5). 
The  steps  are  summarized  as  follows: 

1)  Input  probe  pressure  data,  probe  position,  flow  stagnation 
temperature,  atmospheric  pressure. 

2)  Calculate  ft   and  X  using  the  steps  given  in  Section  A2.U. 

3)  Calculate  K,  with  Eq.  A (22). 

k)     Calculate  K   using  quadratic  interpolation  between  K    values 

-A.  A  « 

1 

given  by  Eq.  A(l9). 

5)   Calculate  X  and  0  using  Eq.  A(23)  and  Eq.  A(2U). 

This  proceedure  was  applied  to  the  calibration  data  as  described  in 
Section  A3. 2  and  listed  in  Table  A2 .   The  results  for  one  survey  are  shown 
for  illustration  in  Table  A7.  At  each  probe  position  in  the  annulus ,  H,  the 
uncorrected  indicated  velocity,  X(UNMOD),  the  measured  actual  velocity, 
X(ACTUAL),  and  the  corrected  velocity,  X(MOD)  are  compared.   The  largest 
error  between  the  corrected  and  actual  values  was  less  than  0.3%  for  all 
points  in  all  surveys  taken. 
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YAW  ANGLE  BALANCE 

FIGURE   Al.      SCHEMATIC    0?   THE    5-H0LE   PROBE 


FIGURE   A2.      VIEW    OF    THE    5-HOLE    PROBE   SHOWING  YAW 

ANGLE   SCALE 
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FIGURE  A3.   SCHEMATIC  OF  THE  APPARATUS  FOR  PROBE 
CALIBRATION  FREE  OF  WALL  EFFECT 
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FIGURE  A4.   PITCH  ANGLE  BLOCK 
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FIGURE  All.   SCHEMATIC  OF  THE  APPARATUS  FOR  PROBE 

CALIBRATION  IN  AN  ANNULUS 


FIGURE  A12.   VIEW  OF  THE  CALIBRATION  APPARATUS 
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Hi-1       %  Hi+1     Hi*2 

FIGURE  Al>.   ILLUSTRATION  OF  QUADRATIC  INTERPOLATION 

METHOD 
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FIGURE  A15.   FLOW  CHART  TO  CORRECT  PROBE  CALIBRATION 

FOR  MEASUREMENTS  IN  AN  ANNULUS 
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APPENDIX  Bi   DETERMINATION  OF  BLOCKAGE 
FACTOR  FOR  FLOW  INTO  THE  ROTOR 
B.l.   INTRODUCTION 

The  analysis  in  Reference   1  established  a 
blockage  factor  for  the  flow  into  the  TRANSX  com- 
pressor rotor  which  was  used  to  correct  inviscid 
mass  flow  calculations.  The  blockage  factor  at 
station  1  was  defined  as, 

Kg  =   m  =  measured  flow  rate   B(l) 
w      inviscid  flow  rate 

and  KB  was  shown  to  be  related  to  the  displacement 
thickness  of  the  boundary  layer.  Preliminary  measure- 
ments  were  made  ahead  of  the  TRANSX  compressor  ro- 
tor to  determine  the  boundary  layer  thickness  at  the 
outer  wall  [Reference   l].  From  the  measurements, 
ignoring  the  hub  wall  boundary  layer,  and  assuming 
that  the  inviscid  flow  field  was  described  correct- 
ly, the  result  obtained  was, 

KB  =  0.989 
Then  the  flow  into  the  rotor  was  described  by  the 
inviscid  flow  calculation  for  a  total  flow  rate 
(w)  which  was  related  to  the  measured  flow  rate 

(w  )  by  the  relation 
m 

*  =  VKBl  • 

The  purpose  of  the  present  study  of  the  flow 
into  the  rotor  was  to  measure  the  inviscid  flow 


94 


field,  as  well  as  the  boundary  layer  region,  and  to 
relate  the  measurements  to  the  analytical  descrip- 
tion. It  was  necessary  therefore  to  change  the  defin- 
ition of  blockage  factor,  from  that  involving  only 
boundary  layers,  to  obtain  the  desired  relationship. 
In  the  following  sections,  necessary  parts  of  the 
analysis  given  in  Reference   1  are  repeated,  and  the 
method  is  given  for  relating  the  analytical  distribu- 
tion of  mass  flux  and  flow  angle  to  the  distributions 
measured  using  the  5-hole  probe. 
B.2.   ANALYTICAL  BACKGROUND  AND  METHOD 

Referring  to  Fig.  Bl,  and  assuming  that  stream- 
lines are  axisymmetric,  continuity  at  station  1 

gives 

r 
w  =  /   °/VcosX«  (2iTr)  dr  B(2) 

ri 

Non-diraensionalizing  Eq.  B(2), 

r2  r   2 

W  -      w       =  r  ro    f V    cos\  d(£)    B(3) 

*  "  *tVro        rf   f*\ 

The  definition  of  the  total  flow  function  is  given 
in  Reference  7   as        , 


$  =  -^r-   =  X(  1-X2)  Y~1,  where 
ft  t 

X  =  V/V.  and  V.  =  V2c  T.  . 
t      t      p  t 

Writing  R  =  r/  ,  Eq.  B(3)  taken  between  R,  and  R2 

o 
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becomes 

4 

W  =  S        I  cos\  d(R2)  • 

4 

Introducing  the  values  of  the  flow  function  and  flow 
angle  at  the  outer  wall,  indicated  by  a  Ho"  subscript, 
the  above  equation  becomes 

.2 


R2 


|     cosX  A(~2 


W  =  £cosX0/  2   -|-   -^  d(R^) 
Rx      °         ° 

In  the  TRANSX  compressor  X      -   0,  therefore 

2 
R. 


W  -  *0  S   2  (-§-)  cosX  d(R2) 
Rl 

The  assumption  is  then  made,  in  Reference  1, 
that  the  pattern  of  streamlines  does  not  depend  on 
the  flow  rate,  that  1|  =  j^/i^is  a  function  of  R  only 
and  cos\  is  another  function  of  R.  Based  upon  func- 
tions derived  from  a  computer  solution  of  the  com- 
pressible flow,  for  R,  =  R.  and  R0  =  R  ,  it  is  shown 

e  1     1        2     0 

that 

2 


m) 


f 


R°-         2 

=  S        I  cosX  d(R^)  =  0.70836  B(ifa) 
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where,  using  Eq.  B(4), 

fft  =  $0^  .  B(4b) 

Therefore  the  "theoretical**  flow  rate  is  proportional 
to  the  value  of  the  analytical  flow  function  at  the 
outer  wall. 

In  order  to  relate  the  analytical  description  of 
the  flow  to  the  measured  distributions  of  the  flow 
function  and  angle,  the  "theoretica-1"  mass  flux  was 
integrated  and  set  equal  to  the  integral  of  the 
measured  distribution  of  the  mass  flux  over  that  re- 
gion of  the  flow  where  the  measurements  were  valid. 
The  probe  measurements  were  considered  to  be  free  of 
"immersion"  or  "wall"  effects  from  20$  to  80$  of  the 
hub-to-tip  survey.  Therefore,  using  Eq.  B(*f) 

2  2 

Rp  Rp 

/  «      cosX  d(R2)    =  ffi     /  j  cosX  d(R2)        B(5) 

r\  ITi  111  O  r\ 

Rl  Rl 

defines  a  value  of  the  analytical  flow  function  at  the 
outer  wall  which  corresponds  to  the  measured  flow 
conditions.  Equation  B(^-b)  then  gave  the  "theoreti- 
cal" flow  rate  corresponding  to  the  measured  flow  rate, 
and  Eq.  B(l)  was  again  used  to  evaluate  a  blockage 
factor. 
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B.3.   PROCEDURE 

The  analytical  distribution  of  the  total  flow 
function  [Ref.   l]  is  represented  as 

(1-R2)  =  a1(l-f)  Vmjd-i),  B(6) 


where  a,,  n,,  m,  are  coefficients,  listed  in  Table  Bl, 
which  were  evaluated  in  Reference   1  from  a  computer 
solution  using  a  curve  fitting  technique,  f  is  un- 
known for  a  given  R  and  must  be  determined  as  follows: 

1)  Choose  an  initial  value  for  $  which  does 
not  exceed  unity.  Take 

Jx   =  1-.03(1-R^)2,  R1   given. 

2)  Set  Y  =  1-Ii  ,  i=l,  2 n 

o       nn 

3)  Calculate  Rt  =  1-a^  -m-jY 

2 

*0   Test  for  convergence  of  calculated  R  ; 


2    2 
R^  -  R7 


<£  ,  £=  convergence  limit. 
If  yes,  the  solution  is  complete. 
If  no,  continue  to  step  5» 
5)   Calculate  the  slope 

2  "I"1 

d(R.)/dy  =  -a,n.,Y     -m,  ,  and  calculate 


Ii+1  -  li  - 


"2   2 
1   l 


d(Rf)/dy 


6)   Repeat  steps  (2)  through  (5)  until  conver- 
gence occurs. 
Using  this  procedure  the  value  of  |>  could  be  calculated 
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for  each' probe  position  R. 

A  quantity  K  was  defined  such  that 

K  =  /     $  cos\.d(R2)  .  B(7) 

2 

Ri 

K  was  determined  for  a  given  R,  and  R2  using  numer- 
ical integration  applied  over  unequal  intervals. 
R,  was  chosen  to  be  the  first  point  greater  than  20$ 
of  the  distance  from  hub  to  tip  and  R2  was  chosen 
to  be  the  last  point  before  80$  of  the  distance  from 
hub  to  tip. 

The  analytical  distribution  for  the  flow  angle 
(A.)  used  in  Eq,  E(^-a)  was  evaluated  in  Reference   1 

as 

n2 
cosX  =  a2  (1-R2)   +m2  (1-R2)  B(8) 

where  the  coefficients  are  given  in  Table  Bl. 

Using  Eq.  B(5),  Eq.  B(^b)  and  Eq.  B(3)  with 
the  perfect  gas  equation  of  state  in  Eq.  B(7), 

f  KB?  roVt   Ro 

Ko  =  0?k    ?         Pt  I.  oosXm  d(R  >       B(9) 

m  G  t      D£ 

Rl 

p. i  X      and  5  are  given  from  measurements  at  each 
rz       m     *m 

point  in  the  probe  survey  and  all  other  quantities  in 
Eq.  B(9)  are  constant  for  the  survey.  Since  slight 
variations  in  the  flow  rate  were  noted  during  the  sur- 
vey, a  correction  was  introduced.  The  measured  flow 
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rate.  w  ,  in  Eq.  B(9)  was  the  average  of  data  reduced 

m 

from  scans  of  fixed  instrumentation  taken  before  and 
after  each  survey.  The  value  of  w  corresponded  to 
a  differential  pressure,  Ah  ,  across  a  flow  nozzle. 
Since  the  flow  rate  at  any  time  was  proportional  to 
the  square  root  of  this  differential  pressure,  Ah, 
the  value  of  Ah  was  recorded  at  each  survey  point. 
Rearranging  Eq,  B(9),  with  the  correction  introduced 
into  the  integration,  the  blockage  factor  is  given 
by 


K 


K  w  R_T . 

-   o  m  G  t 

?  wroVt 


R^  /Ah  ' 


A 


S   2   Vt-u^  P+§  cosX  d(Rfc) 
o  Ah  Kt*m    m 


nl 


B(10) 


Eq,  B(6),  B(7)  and  B(10)  were  programmed  to  calcu- 
late the  blockage  factor  for  flow  at  the  rotor  inlet 
from  probe  survey  measurements. 

The  correction  for  slight  variations  in  the  flow 
rate  was  also  included  in  comparing  the  measured  and 
analytical  distributions  of  the  mass  flux  (or  flow 
function).  The  left  hand  side  of  Eq.  B(5)  is  written, 
with  the  addition  of  the  flow  correction,  as 
.2 

B(ll) 


R! 


K      =     S        § 

m  "    P  xm 


l 


Ahr  2 

Air  cosXm  d(R  ) 


so   that  with  Eq.    B(7).    Eq.    B(5)    becomes 

K      =   $   K    . 
m        *o   o 


100 


The  analytical  flow  function,  which  was  included 

on  plots  of  the  measured  flow  function,  was  obtained 

from 

I(R)  =  I  (R)  (j^J  B(12) 

o 

where  ¥  was  obtained  from  Eq.  B(6),  K  was  obtained 

*■  m 

from  Eq.  B(ll)  and  K  was  already  known  from  Eq.  B(7). 
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TABLE  Bl.   CURVE  FIT  COEFFICIENTS 
FOR  EQ.  B(6)  AND  EQ.  B(8) 


EQUATION 

i 

a. 

1 

n. 

1 

m. 

1 

B(6) 
B(8) 

1 
2 

1.79864 
0.98927 

0.38582 
3.81679 

22.3048 
0.47149 
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